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ABSTRACT 

Long-term records of sunspot number and concentrations of cosmogenic ra¬ 
dionuclides (lOBe and 14G) on the Earth reveal the variation of the Sun’s mag¬ 
netic activity over hundreds and thousands of years. We identify several clear 
periods in sunspot, lOBe, and 14G data as 1000, 500, 350, 200 and 100 years. We 
found that the periods of the hrst five spherical harmonics of the slow magnetic 
Rossby mode in the presence of a steady toroidal magnetic field of 1200-1300 G 
in the lower tachocline are in perfect agreement with the time scales of observed 
variations. The steady toroidal magnetic field can be generated in the lower 
tachocline either due to the steady dynamo magnetic field for low magnetic dif- 
fusivity or due to the action of the latitudinal differential rotation on the weak 
poloidal primordial magnetic field, which penetrates from the radiative interior. 

The slow magnetic Rossby waves lead to variations of the steady toroidal mag¬ 
netic field in the lower tachocline, which modulate the dynamo magnetic held 
and consequently the solar cycle strength. This result constitutes a key point 
for long-term prediction of the cycle strength. According to our model, the next 
deep minimum in solar activity is expected during the hrst half of this century. 
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1. Introduction 


Solar activity, which is manifested through huge energy releases in terms of coronal 
mass ejections and solar flares, basically determines the plasma conditions o f near-Earth 
and interplanetary space, and possibly influences Earth’s climate flHaighl 120071) . Therefore, 
the prediction of solar activity is of paramount importance for our future scientihc- technical 
strat egy. The activity undergoes 11-year oscillations kno wn as Schwabe cycl es (jSchwabe 


18441 ) generally interpreted in terms of dynamo moc 
long-term records of sunspot number flGleissberdll939l: 


mogenic radionuclides on the Earth fIStuiver and Bray.iuna.s 


els ( Charbonneaul 120051) . However, 


Hill 


977) and concentrations of cos- 


19891: ISuessI Il980l: ISolanki et ah 


20041 : lUsoskin et al.ll2004 ISteinhilber et al.ll2012f) reveal the variation of the Sun’s magnetic 
activity over hundreds and thousands of years. 


The observed long-term variation in solar activity is still a mystery flHathawavi 120101) . 
Three different explanations for the variation have been proposed. The hrst one suggests 
that the nonlinear solar dynamo permits a transition from periodi c oscillations to c haotic 
behavior with a long period modulation of magnetic field strength flWeiss et al.l 119841) . The 
second mechanism considers random fluctuations in the dynamo p arameters, which may 
mimi c the main features of the long-term variations in the activity fIChoudhuri and Karak 
20121) . The third proposal is that the tidal influence of solar system planets may modulate the 
temperature gradie nt and consequently the magnetic storage capacity in the solar tachocline 
flAbreu et al.N2012f) . The hrst two mechanisms may account for the occurrence of Grand 
Minima, but fail to reproduce the observed long-term modulation of solar cycles. The third 
one, the planetary hypothesis (a long sta nding problern), ha s been intensively discussed in 
the last years by the scientihc commun ity flGharbonneaull2013[). but has been criticized due to 
the in appropriate stati stical analysis flGameron and Schiisslerl 120131: iPoluianov and U.soskin 
20141 ) or data artefact (IGauquoin et al.l 120141 ). 


The long-term variation of solar activity is probably connected to the tachocline, a thin 
layer below the solar convection zone, which is a very irnporta nt place for the solar dynamo 
and angular momentum redistribution ( Spiegel and Zahnlll992 ). The tachocline is bel ieved to 
consi st of a lower strongly stratihed and an upper overshoot weakly stratihed layers flGilman 
20001) . The convective motions from the convection zone penetrate the upper overshoot layer, 
while the lower layer is m uch more stable and the ordinary magnetohydrodynamic (MHD) 
shallow water equations of iGilmaru (120001) can be safely applied there. 


The conservation of total vorticity in hydrodynamic shallow water equations leads 
to the appearance of Rossby or planetary waves in the rotating frame, which govern the 
large-scale dynamics of the Earth’s atmosphere (and oceans) causing the formation of cy¬ 
clones/anticyclones and hurricanes. Inclusion of horizontal magnetic field splits the ordinary 
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hydrody namic Rossby waves into fast and slow magnetic Rossby waves in a magnetized 
plasma fjZ agar ashy ili et al.l 120071. 1200911 . Fast magnetic Rossby waves are responsible for 
the observed Rieger type-periodicity (155-160 da ys') in solar activity after the interaction 
with the dynamo magnetic held in the tachocline fjZaqarashyili et al.l 1201011 . However, slow 
magnetic Rossby wayes conld have very long time scales compared to the solar rotation and 
solar cycle periodicity. In this letter, we consider the slow magnetic Rossby modes in the 
lower strongly stable layer of the solar tachocline and show their connection with long-term 
variations in solar activity. 


2. Slow magnetic Rossby waves and long periods 


J.Xi OWJ.CI.1. K J. V J. 


We consider the spherical coordinate system (r, 6*,0), where r is the radial coordinate, 
9 is the co-latitnde and 0 is the longitnde. The nnpertnrbed magnetic held is assnmed to 
be toroidal i.e. B = B^{6)e^. Then, the consideration of rigid rotation and a homogeneons 
magnetic held, = Bq = const, l eads to the following d ispersion relation for magnetic 
Rossby wayes in the rotating frame fjZaqarashyili et al.ll2007l) 


+ 


2flom B^m? 2 — n{n + 1) 

- UJ - 

n{n -I- 1) n{n -\- 1) 


= 0 , 


( 1 ) 


where ca is the freqnency of the spherical harmonic, flo is the eqnatorial angnlar velocity, 
R is the distance from the solar centre, p is the medinm density, m and n are toroidal and 
poloidal wavennmbers, respectively. If Bq/{AtipVLqR^) 1, which is easily satished in the 
solar tachocline, then the low order spherical harmonics of slow magnetic Rossby modes are 
described by the dispersion relation 


(jJr, 


= —mflr 




2 — n{n + 1) 


A'npVL^R? 


( 2 ) 


where Umn is the freqnency of the spherical harmonic. Each spherical harmonic with different 
m and n has a different freqnency (and hence period) depending on the plasma parameters; 
mostly on the magnetic field strength, Bq. 

In order to find the correlation between slow magnetic Rossby waves and the long-term 
variation in solar activity, the periods of spherical harmonics with different magnetic field 
strength shonld be compared to the time scales fonnd in observational data. We performed a 
spectral analysis of yearly mean total snnspot nnmber from 1700 to 2013, and fonnd strong 
peaks at 100 and 180 years in the oscillation spectrn m (Fignre 1, npper panel). The pe riods 
can be identified with the well known Gleissberg flGleissberd Il939l: iHathawayi 120101) and 
Sness f)Snesslll980l) cycles, respectively. Longer oscillatory periods cannot be safely defined 
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due to the short length of the data series (~ 310 year s'). To hnd longer periods, we used 


solar activity reconstructed from cosmogenic isotopes fjVonmoos et ah 


2006: Usoskin et ah 


20071) lOBe and 14C measured in the Greenland GRIP ice core and in tre e rings, respectively, 


during the past 10 mill e nnia. Both data sets sh ow different periodicities flUsoskin et al.ll2004J: 


Steinhilber et al.l 120121: iHanslmeier et al.l 120131) . In order to select reliable periods, we put 


forward a novel approach and overlapped the spectra of both data sets (Figure 1, lower 
panel). Then we selected four well defined periods of 1000, 500, 350 and 200 years in both 
data sets. lOBe data also shows power at the periods of 700-800 years, but there are no 
corresponding periods in 14G data. Gombining the power spectra of sunspot numbers and 
cosmogenic isotopes, we dehne the reliable periods in solar activity as ~ 1000, ~ 500, ~ 350, 
~ 200 and ~ 100 years. 


Using the plasma parameters at the tachocline (Gq =2.7T0“® s“^, p =0.2 g-cm“^, R = 
5T0^° cm) and the magnetic field strength of 1200 G, we find that the periods of the modes 
m=l, n =2, 3, 4, 5, 6 as calculated from Eq. (2) are ~ 1170, ~ 470, ~ 260, ~ 170 and 
~ 110 years, respectively, which show perfect agreement between theoretical and observed 
periods (Figure 2). Therefore, we suggest that the observed periods in solar activity are the 
consequence of the first five spherical harmonics of slow magnetic Rossby modes in the lower 
tachocline. 


The magnetic held strength of 1200 G, considered above, is much smaller than the esti¬ 
mated value of dynamo generated magnetic held in the tachocline (10^10^ G). Neyertheless, 
the toroidal component of the reversing dynamo held is zero when averaged over a time 
interval longer than the solar cycle length, because of periodic sign reversals in consecutive 
cycles. Therefore, it probably cannot inhuence the long-period dynamics of slow magnetic 
Rossby waves. However, there are two diherent mechanisms, which may be responsible for 
the appearance of weak steady (n onreversing) toroidal magnetic held in the lower tachocline 
or upper radiative interior. First, iDikpati et al.l (120061) showed that a steady (nonreyersing) 
dynamo for a low enough magnetic dihusivity < 10^ cm^ s“^ may generate steady toroidal 
magnetic held with the strength of ~ 1 kG in the lower tachocline. This is consistent to 
our requirements. Second, the primordial magnetic held, which remained in the radiatiye 
interior after the Sun’s collapse, may penetrate into the tachocline and ahect Rossby waves. 
Both, nonreversing dynamo and the primordial helds m ay be responsible for the observed 
Gnevyshev-Ohl 22-year rule flGnevyshev and Ohll 119481) . which shows that odd sunspot cy¬ 
cles are genera lly stronger than even cycles, due to the break of polarity symmetry of the 
dynamo state flBoyer and Leyvl 119851) : the reyersing dynamo magnetic held is aligned and 
anti aligned with respect to the toroidal component of steady held in consecutive cycles 
leading to the consequence of strong and weak cycles, respectively. Similarly, the periodical 
variation caused by slow magnetic Rossby wayes in the steady magnetic held can inhuence 
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the strength of the dynamo magnetic held and conseqnently the cycle strength. 


The primordial magnetic held strength was estimated as 1000 G deep in the radiative 
interior and 1 G jnst beneath the tachocline flGough and Mclntvrelll9981) . This is the poloidal 
magnetic held, which probably has negligible inhnence on the cycle strength. Helioseismology 
shows that the radiative interior rotates almost nniformly, therefore the magnetic held will 
remain poloidal in most of its part. On the other hand, the convection zone and the tachocline 
have diherential rotation expressed as 

G = r2o(l — S 2 cos^ 0 — S 4 cos^ 6 *), (3) 


rose valnes at the solar snrface 


where S 2 and S 4 are parameters determined by o bservations w , 

(and in the middle tachocline) are abont 0.14. ISchou et al.l fjl998l) constrncted (see Figure 


7 in that paper) the radial dependence of solar rotation at diherent latitudes along the 
whole convection zone and radiative interior as inferred by helioseismology. The authors 
determined the center of the tachocline at 0.7-0.71 Rq with the negligible thickness of 0.05 
Rq, therefore the lower bo undary of th e tach ocline is located somewhere at 0.68 Rq. A 
careful look at Figure 7 of ISchou et al.l (1199811 shows that G is 445 nHz at 0=0 and 415 
nHz at 0=60 degree near the tachocline center (0.7 Rq), which gives the differential rotation 
parameter S2=0.1. This is close to the main value of differential rotation in the tachocline. 
On the other hand, G is 440 nHz at 0=0 and 435 nHz at 0=60 degree near the lower boundary 
of the tachocline (at 0.68 Rq), which gives §2=0.015. Therefore, a small differential rotation 
is still present in the lower tachocline. This differential rotation may eventually lead to 
the generation of a toroidal component from the poloidal primordial held, which penetrates 
from the upper part of the radiative envelope into the lower part of the tachocline. The 
induction equation allows to estimate the generated toroidal held as B^/Iq ^ Bg{dfl/d0), 
where to is the characteristic time (taken as the Alfven time R^/Anp/Bg) and Bg is the 
poloidal component of the primordial magnetic held. Then the toroidal magnetic held is 


Bff, ~ Bq cos 0 sin 0, (4) 

where Bq is 4300 G for § 2 = 0.01 and 1500 G for §2=0.0035. The strength of the estimated 
magnetic held coincides with the value required for slow magnetic Rossby waves to explain 
the observed periodicity. It should be noted that the equipartition value of steady toroidal 
magnetic held with the kinetic energy of turbulent pumping from the convection zone ex¬ 
pressed as Beq ~ V\/Anp, where v is the turbulent velocity of convective cells (0.01 km-s“^ 
near the base of the convection zone according to the standard mixing length theory of 
convection), is also estimated as 1500 G. 


Therefore, the strength of the steady toroidal magnetic held in the lower tachocline 
generated either by steady dynamo for low magnetic dihusivity or by diherential rotation 
acting on the primordial poloidal component can be safely taken as 1000-1500 G. 
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The solar cycle strength is dehned by magnetic energy of dynamo layer, therefore it is 
connected to the square of magnetic held strength. Longitudinally averaged magnetic energy 
of toroidal magnetic held is proportional to {{Bd + + {^BdBs)^ + (B^)^, 

where brackets denote averaging over (j) from 0 to 27r, B^ is the dynamo generated reversing 
magnetic held and Bg is the steady nonreversing (primordial or steady dynamo) magnetic 
held. Spherical harmonics of slow magnetic Rossby waves lead to the toroidal dependence of 
steady magnetic held corresponding to the wave number, m. For example, m = 1 mode leads 
to the dependence of cos0 in the steady held. Then the longitudinally averaged magnetic 
energy is proportional to Bj + B‘^q/2, where Bgo is the amplitude of steady held. Bsq has 
temporal variations due to the slow magnetic Rossby waves as discussed above, therefore 
it can be responsible for the long period modulation of longitudinally averaged magnetic 
energy and hence the solar cycle strength. 


The latitudinal diherential rotation in the lower part of the tachocline (or in the upper 
part of the radiative interior), even if small, may l ead to magnetohydrodynamic instabilitie s 
in joint action with the toroidal magnetic held flGilman and FoxI 119971: ICallv et al.l l2003l) . 
In order to hnd unstable ha.rmoni cs, we use the general technique of Legendre polynomial 
expansion flLonguet-HigginsI 119681) . 2D shallow water MHD equations in the presence of 
diherential rotatio n (Eg. |3l) and a toroid al magnetic held (Eq. H]) can be written in the 
rotating frame as fjZaqarashvili et al.ll2010l) 


(Drf - u)L^ + (2 - - /i')])^ - - 6/i/32<F = 0, (5) 


where 


(Drf -a;)<F = /iT, 


r ^ n . 2 \ 9 

l-fj? 


( 6 ) 


is the Legendre operator, T and $ are the stream functions for velocity and magnetic held 
perturbations, /x = cos 6 and 


.2 „ ,.4 ^ q2 ^0 


Grf(/i) S2/i" S4/, (n 

We expand T and $ in inhnite series of associated Legendre polynomials 

00 00 

n=m n=m 

substitute them into Eqs. ([5])-([6]) and, usi ng a recurrence relation o f Legendre polynomials, 
obtain algebraic equations as inhnite series flZaqarashvili et al.ll2010l) . The dispersion relation 
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for the infinite number of harmonics can be obtained when the inhnite determinant of the 
system is set to zero. In order to solve the determinant, we truncate the series at n = 75 and 
solve the resulting polynomial in oj numerically. The frequencies of different harmonics can 
be real or complex giving the stable or unstable character of a particular harmonic. Then 
differential rotation parameters S 2 = S 4 = 0.01 (as estimated in the previous paragraph) and 
the magnetic field strength Bo = 1200-1300 G lead to periods of unstable harmonics - which 
are close to the observationally obtained and analytically predicted values (Figure 2). 


3. Discussion and Conclusions 


It is shown in this letter that the spherical harmonics of slow magnetic Rossby waves 
have very long periods in the presence of a steady toroidal magnetic field near the lower 
boundary of the tachocline. The steady toroidal field can be generated in the lower tachocline 
either due to a steady dynamo in the case of low magnetic diffusivity (IDikpati et al.l 120061) 
or due to action of latitudinal differential rotation on poloidal primordial magnetic field. 
The strength of steady dynamo magnetic hel d in the lower tacho cline is ~ 1 kG in the case 
of low magnetic diffusivity (< 10^ cm^ s“^) flDikpati et al.ll2006l) . The poloidal primordial 
magnetic held may penetrate the lower part of the tachocline from the radiative interior 
and generate the toroidal component due to the small latitudinal diheren tial rotation, whic h 
is present near the interface of the radiative interior and the tachocline fjSchou et al.lll998h . 
The diherential rotation rate of the order of S 2 = 0.003 — 0.01 leads to the equipartition value 
of toroidal magnetic held of 1500-4000 G. In both cases of steady dynamo and primordial 
magnetic helds, the periods of spherical harmonics of slow magnetic Rossby waves correspond 
to the observed long periods in sunspot number and radionuclide data. We believe that 
this coincidence has reliable physical ground: slow magnetic Rossby modes cause periodical 
variations of steady nonreversing toroidal magnetic held in the lower tachocline, which in 
turn modify the dynamo magnetic held and hence the solar cycle strength. The lower part 
of the tachocline is stable contrary to the upper overshoot layer as the convection does not 
penetrate there. Therefore, the steady toroidal magnetic held and large-scale slow magnetic 
Rossby waves may survive over long time without signihcant inhuence. 






















Rossby waves along the toridal direction is 


Vph 


Bl 2 — n{n + 1) 
47rpf2o-R 2 


( 8 ) 


which for n = 2 and i?o=1200 G has the value of 8 cm s“^. The speed is much slower 
compared to the surface speed of solar rotation, therefore it is hard to catch the migration 
by observations. 


Slow magnetic Rossby waves, which reproduce the observed long periods in solar activity, 
can be used for the long-term forecasting of the activity level. We £t the sum of hve sinusoidal 
functions (with periods of 1000, 500, 350, 200 and 100 years), which are actually the spherical 
harmonics of slow magnetic Rossby modes, to the data of yearly mean total sunspot number 
from 1700 to 2013 (Figure 3). The resulting curve hts rather well the long-term behavior of 
the activity. Our model predicts the next deep minimum of solar activity around year 2030, 
which means that the next two solar cycles will be rather weak, but then the activity will 
start to increase again. 


We also fit the sum of four sinusoidal functions (1000, 500, 350, 200 years) to the 100- 
year averaged sunspot number for the last three millennia, reconstructed from 14C data 
(Figure 4). The £t is generally good except near the Sporer minimum, where a remarkable 
anti phase relation is seen. This disagree ment is probab l y cau sed due to two distinct modes 
in solar activity as recently suggested by lUsoskin et al.l (120141) . 


The theory we have described provides a simple view of the coupling between the large- 
scale dynamics of the radiative interior and the tachocline. Dynamic interaction between 
rotation and steady (steady dynamo or primordial) toroidal magnetic held in the upper 
part of the radiative interior and the tachocline inevitably results in long periodicities due 
to slow magnetic Rossby modes, which accordingly modify the reversing dynamo magnetic 
held causing the observed long-term variations in solar activity. The scenario is a starting 
point for future sophisticated models with comprehensive study (including intense numerical 
simulations), which will enhance the accuracy of long-term solar activity forecast. The 
theory may also provide a hint for observed long-period modulations of magnetic activity in 
solar-like stars. 
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Period, years 



Fig. 1.— Upper panel: Periodogram of yearly mean snnspot nnmbers. Large power valnes 
are concentrated aronnd 180 and 100 years. Power is given in arbitrary nnits. Yearly 
mean snnspot data is downloaded from the World Data Center-SILSO, Royal Observatory 
of Belginm, Brnssels. Low er panel: periodogram s of reconstrncted sol ar activity based on 
cosmogenic isotopes lOBe (IVonmoos et al.ll2006l ) (black line) and 14C flU.so.skin et aLl 1200711 
(red line). Large power valnes in both data are concentrated aronnd 1000, 500, 350 and 200 
years. Power is given in arbitrary nnits. 
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Fig. 2.— Periods of solar activity and slow magnetic Rossby modes with different poloidal 
wave number, n. Blue and Red circles indicate the periods found by periodogram analysis of 
sunspot (Figure 1, upper panel) and radionuclide (Figure 1, lower panel) data, respectively. 
Blue and magenta asterisks indicate the periods of slow magnetic Rossby wave harmonics 
calculated from Eq. (2) for m=l and n=2, 3, 4, 5, 6 for magnetic held strength of 1200 
and 1300 G, respectively. Blue and magenta triangles indicate the periods of the hrst hve 
unstable harmonics calculated by the Legendre polynomial expansion for magnetic held 
strength of 1200 and 1300 G, respectively. Here m=l and the diherential rotation parameters 
are S 2 = S 4 = 0 . 01 . 
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1700 1750 1800 1850 1900 1950 2000 2050 
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Fig. 3.— Forecast of future solar cycle strength using magnetic Rossby wave theory. Red 
and blue lines show yearly and 10-year averaged sunspot numbers, respectively. Magenta 
line shows the £t of the sum of 5 sinusoidal functions with periods of 1000, 500, 350, 200 
and 100 years. 
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Fig. 4.— Fit of magnetic Rossby wave harmonics to reconstructed sunspot number for the 
last three millennia. Red line shows 100-year averaged sunspot number reconstructed from 
14C data fjUsoskin et al.ll2014l ). Purple line shows the £t of the sum of 4 sinusoidal functions 
with periods of 1000, 500, 350 and 200 years, which are the hrst four harmonics of slow 
magnetic Rossby waves. 






